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Abstract 

We consider principal pivot transform {pivot) on graphs. We define a natural 
variant of this operation, called dual pivot, and show that both the kernel and 
the set of maximally applicable pivots of a graph are invariant under this oper- 
ation. The result is motivated by and applicable to the theory of gene assembly 
in ciliates. 

Keywords: principal pivot transform, algebraic graph theory, overlap graph, 
gene assembly in ciliates 



1. Introduction 

The pivot operation, due to Tucker partially (component-wise) inverts 
a given matrix. It appears naturally in many areas including mathematical 
programming and numerical analysis, see [l7| for a survey. Over F2 (which is 
the natural setting to consider for graphs), the pivot operation has, in addition 
to matrix and graph interpretations also an interpretation in terms of delta 
matroids [1]. 

In this paper we define the dual pivot, which has an identical effect on graphs 
as the (regular) pivot, however the condition for it to be applicable differs. The 
main result of the paper is that any two graphs in the same orbit under dual 
pivot have the same family of maximal pivots (cf. Theorem [TO]) . i.e., the same 
family of maximally partial inverses of that matrix. This result is obtained by 
combining each of the aforementioned interpretations of pivot. 

This research is motivated by the theory of gene assembly in ciliates Q, 
which is recalled in Section [T] Without the context of gene assembly this main 
result (Theorem I16p is surprising; it is not found in the extensive literature 
on pivots. It fits however with the intuition and results from the string based 
model of gene assembly and in this paper we formulate it for the more 
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general graph based model. It is understood and proven here using completely 
different techniques, algebraical rather than combinatorial. 



2. Notation and Terminology 

The field with two elements is denoted by IF2- Our matrix computations will 
be over F2. Hence addition is equal to the logical exclusive-or, also dcinoted 
by ®, and multiplication is equal to the logical conjunction, also denoted by 
A. These operations carry over to sets, e.g., for sets A,BCV and x G V, 
X G A ® B iS {x & A) ® {x G B). 

A set system is a tuple M = {V, D), where V is a finite set and D C 2^ is 
a set of subsets of V. Let min(I?) (max(£)), resp.) be the family of minimal 
(maximal, resp.) sets in D w.r.t. set inclusion, and let min(M) = (y, min(Z))) 
(max(M) = (y, max(Z))), resp.) be the corresponding set systems. 

Let y be a finite set, and A be a ^ x F-matrix (over an arbitrary field), i.e., 
A is a matrix where the rows and columns of A are identified by elements of V. 



Therefore, e.g., the following matrices with V = {p,q} are equal: 

Q 



1 1 
1 



and ' 2^ ][ J • X CV, the principal submatrix of A w.r.t. X is denoted 

by ^[-'i^], i.e., A[X] is the X x X-matrix obtained from A by restricting to rows 
and columns in X. Similarly, we define A\X = A[V\X]. Notions such as matrix 
inversion A~^ and determinant det(A) are well defined for V x F-matrices. By 
convention, det(A[0]) = 1. 

A set X C V is called dependent in A iff the columns of A corresponding 
to X arc linearly dependent. We define Va =■ {I, D) to be the partition of 2^ 
such that D (/, respectively) contains the dependent (independent, respectively) 
subsets of V in A. By convention, g /. The sets in max(/) are called the 
bases of A. 

We have that Va = (I, D) is uniquely determined by max(/) (and the set 
V). Similarly, Va is imiqucly determined by min(_D) (and the set V). These 
properties are specifically used in matroid theory, where a matroid may be 
described by its independent sets (V,/), by its family of bases (V,max(/)). or 
by its circuits (V,min(L')). Moreover, for each basis X G niax(/), \X\ is equal 
to the rank r of A. 

We consider undirected graphs without parallel edges, however we do allow 
loops. For a graph G — (V, E) we use V{G) and E{G) to denote its set of 
vertices V and set of edges E, respectively, where for a; G V", {x} G ii^ iff x has 
a loop. For X C V", we denote the subgraph of G induced by X as G'[A]. 

With a graph G one associates its adjacency matrix A{G), which is a x V- 
matrix (a„,„) over F2 with ciu.v 

= 1 iff {u, v} G E. The matrices corresponding to 
graphs are precisely the symmetric F2-matrices; loops corresponding to diagonal 
I's. Note that for X<ZV, A{G[X]) = {A{G))[X]. 



g 

1 



2 



Over F2, vectors indexed by V can be identified with subsets of V, and a 
V X TZ-matrix defines a linear transformation on subsets of V. The kernel (also 
called null space) of a matrix A, denoted by ker(A) is determined by those linear 
combinations of column vectors of A that sum up to the zero vector 0. Working 
in F2, we regard the elements of ker{A) as subsets of V. Moreover, the kernel 
of A is the eigenspace Eq{A) on value 0, and similar as ker(A), the elements 
of the (only other) eigenspace Ei (A) ^ {v ^ V \ Av — v} on value 1 are also 
considered as sets. 

We will often identify a graph with its adjacency matrix, so, e.g., by the de- 
terminant of graph G, denoted by det G, we will mean the determinant det A{G) 
of its adjacency matrix computed over F2. In the same vein we will often simply 
write ker(G), Ei{G), Vg, etc. 

Let Vg = (I,D) for some graph G. As G is a V x V^-matrix over F2, 
we have that X e D iS there is a 5* C X with S £ ker(G)\{0}. Moreover, 
min(D) — min(ker(G)\{0}) and ker(G) is the closure of min(Z)) under © (i.e., 
min(Z?) spans ker(G)). Consequently, min(_D) uniquely determines ker(G) and 
vice versa. As min(_D) in turn uniquely determines Vg, the following holds. 

Corollary 1. For graphs Gi and G2, ker(Gi) = ker(G2) iff the families of 
bases of Gi and of G2 are equal. 



3. Pivots 



In general the pivot op eration can be studied for matrices over arbitrary 
fields, e.g., as done in [l7|. In this paper we restrict ourselves to symmetric 
matrices over F2, which leads to a number of additional viewpoints to the same 
operation, and for each of them an equivalent definition for pivoting. Each of 
these definitions is known, but (to our best knowledge) they were not before 
collected in one text. 



Matrices. Let A be a x T^-matrix (over an arbitrary field) , and let X CV be 
such that A[X] is nonsingular, i.e., det ^ 0. The pivot of A on X, denoted 



hy A* X , is defined as follows, see \li 
Then 



Let A = 



p 


Q 


R 


s 



with P ^ A[X]. 





-P-'Q 


RP-^ 


S-RP-'Q 



Matrix {A * X) \X = S — RP~^Q is called the Schur complement of X in A. 

The pivot is sometimes considered a partial inverse, as A and A* X are 
related by the following characteristic equality, where the vectors xi and yi 
correspond to the elements of X. In fact, this formula defines A* X given A 
and X [l7|. 



A 



Xi 
X2 



yi 
y2 



iS A*x 



yi 

X2 



Xi 

2/2 



(1) 



3 



Note that if detA ^ 0, then A^V^ A~'^. By Equation ((J) we see that a 
pivot operation is an involution (i.e., operation of order 2), and more generally, 
if {A * X) * y is defined, then A* {X (BY) is defined and they are equal. 

The following fundamental result on pivots is due to Tucker [3l (see also 0, 
Theorem 4.1.1]). It is used in [sj to study sequences of pivots. 

Proposition 2 ([18]). Let A be a V x V -matrix, and let X C V be such that 
detA[X] ^ 0. Then, forY C V, det{A*X)[Y] = det A[X ® Y] / det A[X] . 

It may be interesting to remark here that Proposition [5] for the case Y = V\X 
is called the Schur determinant formula and was shown already in 1917 by Issai 
Schur, see [lij . 

It is easy to verify from the definition of pivot that A* X is skew-symmetric 
whenever A is. In particular, if G is a graph (i.e., a symmetric matrix over F2), 
then G * X is also a graph. From now on we restrict our attention to graphs. 

Delta Matroids. Consider now a set system M = (V, D). We define, for X C_ V, 
the twist M*X = {V,D*X), where D * X ^ {Y (B X \ Y e D}. 

Let G be a graph and let Mg = {^{G), Dq) be the set system with Dq = 
{X C V{G) I det G[X] — 1}. It is easy to verify that G can be (re)constructed 
given Mg- {u} is a loop in G iff {u} G Dq, and {u,v} is an edge in G iff 
i{u,v} e Dg) © {{{u} e Dg) a {{v} G Dg)), see 0, Property 3.1]. In this way, 
the family of graphs (with set V of vertices) can be considered as a subset of 
the family of set systems (over set V). 

Proposition [5] allows for another (equivalent) definition of pivot over F2 . 
Indeed, over F2, we have by Proposition [H det(^ * X)[Y] = detA[X © Y] for 
all y C y assuming A* X is defined. Therefore, for A4g*x we have Dg*x = 
{Z I det((G * X)[Z]) = 1} = {Z\ dct{G[X ® Z]) ^ 1} ^ {X ®Y \ det(G[r]) = 
1} = Dg * X, see [l|. Hence Mg * X = Mg*x is an alternative definition of 
the pivot operation over F2. 

It turns out that M.g has a special structure, that of a delta matroid, allowing 
a specific exchange of elements between any two sets of DG^ see [if. However, 
not every delta matroid M has a graph representation, i.e., M may not be of 
the form M.g for any graph G (a characterization of such representable delta 
matroids over F2 is given in 0). 

Example 3. Let G be the graph depicted in the upper-left corner of Fig- 



ure □ We have A{G) 



/ 1 1 1 \ 

110 1 

10 1 

V 1 1 1 y 



This corresponds to M.g 



[{p,<l,r,s},DG), where 

Dg = {0, {q}, [P, g}, [P, r}, {p, s}, {g, s}, {r, s}, {p, q, s}, {p, q, r}, {q, r, s}}. 
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Figure 1: The orbit of G under pivot. Only the elementary pivots are showrn. 

For example, {p,q} G Dg since det{G[{p,q]]) — det ^ ^ ^ = 1. Then 

Dg * {P, q] = {0, {p}, {r}, {s}, {p, q}, {p, s}, {q, r}, {q, s}, {p, r, s}, {p, q, r, s}}, 
and the corresponding graph is depicted on the top- right in the same Figure [TJ 
Equivalently, this graph is obtained from G by pivot on {p,q}. Also note that 
we have Dg * {p,s} = Dg, and therefore the pivot of G on {p, s} obtains G 
again. The set inclusion diagrams of M-g ^-nd MG*{p,q} ^-re given in Figured 

Graphs. The pivots G * X where X is a minimal element of M.g\{^} w.r.t. 
inclusion are called elementary. It is noted in that an elementary pivot X 
corresponds to either a loop, X — {u} £ E{G), or to an edge, X = {u,v} £ 
E{G), where both vertices u and v are non-loops. Moreover, each Y G Mg can 
be partitioned Y = XiU- ■ -UX,, such that G*Y = G*{Xi(B- ■ -^Xn) = (• • • (G* 
Xi) • • • * Xn) is a composition of disjoint elementary pivots. Consequently, a 
direct definition of the elementary pivots on graphs G is sufficient to define the 
(general) pivot operation. 

The elementary pivot G*{u} on a loop {u} is called local complementation. It 
is the graph obtained from G by complementing the edges in the neighbourhood 
NgIu) = {v e V \ {u,v} e E{G),u ^ v} oi u in G: for each v,w e Ng{u), 
{v, w} e E{G) iff {v, w} ^ E{G * {u}), and {v} G E{G) iff {v} ^ E{G * {u]) 
(the case v = w). The other edges are left unchanged. 

The elementary pivot G * {u, u} on an edge {u, w} between distinct non- 
loop vertices u and v is called edge complementation. For a vertex x consider 
its closed neighbourhood Nq{x) — Ng{x) U {x}. The edge {m, w} partitions 
the vertices of G connected to u or u into three sets Vi = N'q{u) \ Nq{v), 

= N^{v) \ iV^(M), V3 = N^{u) n N[.{v). Note that u, w G ^3- 

The graph G*{u, v} is constructed by "toggling" all edges between different 
Vi and Vj-. for {x, y} with x e Vi and y <E Vj {i ^ j): {x, y} G E{G) iff {a;, y} ^ 
£'(G'[{ii, w}]), see Figure[2] The remaining edges remain unchanged. Note that, 
as a result of this operation, the neighbours of u and v are interchanged. 

Example 4. The whole orbit of G of Example[3]under pivot is given in FigureHJ 
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Figure 2: Pivoting {u,v} in a graph. Connection {x,y} is toggled iff z g and y (H Vj witli 
i j- Note tliat u and v are connected to all vertices in V3, these edges are omitted in the 
diagram. The operation does not affect edges adjacent to vertices outside the sets Vi, V2, V3, 
nor does it change any of the loops. 

It is obtained by iteratively applying elementary pivots to G. Note that G*{p, q} 
is defined (top-right) but it is not an elementary pivot. 

4. Dual Pivots 

In this section we introduce the dual pivot and show that it has some inter- 
esting properties. 

First note that the next result follows directly from Equation ([T]|. 

Lemma 5. Let A be a V x V -matrix (over some field) and let X Q V with 
A[X] nonsingular. Then the eigenspaces of A and A* X on value 1 are equal, 
i.e., Ei[A) ^Ei{A*X). 

Proof. We have v e Ei{A) iS Av = v iff {A * X)v = v iS v e Ei{A * X). □ 

For a graph G, we denote G -I- / to be the graph having adjacency matrix 
A{G) + I where / is the identity matrix. Thus, G -I- / is obtained from G by 
replacing each loop by a non-loop and vice versa. 

Definition 6. Let G be a graph and let X C F with det((G-|-/)[X]) = 1. The 
dual pivot of G on X, denoted by G*X, is ((G + /) * X) + I. 

Note that the condition det((G + = 1 in the definition of dual pivot 

ensures that the expression ((G -I- /) * X) + / is defined. The dual pivot may be 
considered as the pivot operation conjugated by addition of the identity matrix 
I. As / + / is the null matrix (over F2), we have, similar as for pivot, that dual 
pivot is an involution, and more generally {G*X)*Y, when defined, is equal to 
Gi{X®Y). 

By Lemma [5l we have the following result. 

Lemma 7. Let G be a graph and let X <Z V such that G*X is defined. Then 
ker(G5X) =ker(G). 
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Proof. Note that Ax = OiS {A + I)x = x. Hence, ker(G) = Ei{G + 1). Since 
/ + / is the null matrix (over F2), we have also ker(G + I) = Ei (G). 

Therefore, we have ker(G*X) = ker(((G + I)*X) + I)= Ei({G + I)*X) = 
Ei{G + /), where we used Lemma [5] is the last equality. Finally, Ei{G + /) = 
ker(G) and therefore we obtain ker(G*X) = ker(G). □ 

In particular, for the case X = V , we have that ker((G + I)^^ + I) — ker(G) 
(the inverse is computed over F2) if the left-hand side is defined. 

Remark 8. By Lemma[7]and Corollary [1] we have that the (column) matroids 
associated with G and G*X are equal. Note that here the matroids are obtained 
from the column vectors of the adjacency matrices of G and G*X; this is not to 
be confused with graphic matroids which are obtained from the column vectors 
of the incidence matrices of graphs. 

We call dual pivot GiX elementary if *X is an elementary pivot for G + I. 
Equivalently, they are the dual pivots on X for which there is no non-empty 
Y C X where G*Y is applicable. An elementary dual pivot *{u} is defined on 
a non-loop vertex u, and an elementary dual pivot *{u, v} is defined on an edge 
{u, v} where both u and v have loops. This is the only difference between pivot 
and its dual: both the elementary dual pivot *{u} and the elementary pivot 
*{u} have the same effect on the graph — both "take the complement" of the 
neighbourhood of u. Similarly, the effect of the elementary dual pivot v} 
and the elementary pivot v} is the same, only the condition when they can 
be applied differs. 

Note that the eigenspaces Eo{G) = ker(G) and Ei{G) have a natural inter- 
pretation in graph terminology. For X C T^(G), X £ Eo{G) iff every vertex in 
V{G) is connected to an even number of vertices in X (loops do count). Also, 
X e Ei{G) iff every vertex in V{G)\X is connected to an even number of ver- 
tices in X and every vertex in X is connected to an odd number of vertices in 
X (again loops do count). 

Example 9. Let G' be the graph depicted on the upper- left corner of Figure [31 



We have A{G') 
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r 


s 




/I 
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1 
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1 








1 




1 
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V 1 


1 


1 


1 


/ 



Note that Ei{G') = {0,{p,r,s}} is of 



dimension 1. We can apply an elementary pivot over p on G' . The resulting 
graph G' * {p} is depicted on the upper-right corner of Figure [31 and we have 

p q r s 



A{G' * {p}) 



/I 1 1 1\ 



Note that the elements of Ei{G') are pre- 



1110 
110 

\ 1 y 

cisely the eigenvectors (or eigensets) on 1 for A{G' * {p}), cf. Lemma [5] The 
graphs G' + 1 (which is G in Example ^ and G' * {p} + I are depicted in the 
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Figure 3: Dual pivot of graph G from Example |3] (G is shown in the lower-left corner). 

lower-left and lower-right corner of Figure |3l respectively. By definition of the 
dual pivot we have G*{p} = (G" + I)*{p} = G" * {p} + 1. 

It is a basic fact from linear algebra that elementary row operations retain 
the kernel of matrices. Lemma [7] suggests that the dual pivot may possibly be 
simulated by elementary row operations. We now show that this is indeed the 
case. Over F'2 the elementary row operations are 1) row switching and 2) adding 
one row to another (row multiplication over F2 does not change the matrix) . The 
elementary row operations corresponding to the dual pivot operation are easily 
deduced by restricting to elementary dual pivots. The dual pivot on a non-loop 
vertex u corresponds, in the adjacency matrix, to adding the row corresponding 
to u to each row corresponding to a vertex in the neighbourhood of u. Moreover, 
the dual pivot on edge {m, v} (where both u and v have loops) corresponds to 
1) adding the row corresponding to u to each row corresponding to a vertex in 
the neighbourhood of v except it, 2) adding the row corresponding to v to each 
row corresponding to a vertex in the neighbourhood of u except u, 3) switching 
the rows of u and v. Note that this procedure allows for another, equivalent, 
definition of the regular pivot: add /, apply the corresponding elementary row 
operations, and finally add / again. 

Note that the dual pivot has the property that it transforms a symmetric 
matrix to another symmetric matrix with equal kernel. Applying elementary 
row operations however will in general not obtain symmetric matrices. 

5. MELximal Pivots 

In Section [3] we recalled that the minimal elements of A^g, corresponding to 
elementary pivots, form the building blocks of (general) pivots. In this section 
we show that the set of maximal elements of g , corresponding to "maximal 
pivots", is invariant under dual pivot. 
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Figure 4: Set inclusion diagram of A4g, J^G*{p,q}^ ^^'^ ^Gi{p} ^'^^ G*{p,q}, and G*{p} 
as given in Examples [3] and [9] 



For Mg = {V,Dg), we define J"g = maxiDa). Thus, for X C V{G), 
X G -Fg iff det G[X] = 1 while det G[y] = for every Y D X. 

Example 10. We continue Example [S] Let G be the graph on the lower-left 
corner of Figure [S] Then from the set inclusion diagram of A^g in Figure [4] 
we see that J-q = {{p, q, s}, {p, q, r}, {q, r, s}}. Also we see from the figure that 

^G*{p..q} = {V}- 

Next we recall the Strong Principal Minor Theorem for (quasi-) symmetric 
matrices from — it is stated here for graphs (i.e., symmetric matrices over 

F2).E1 

Proposition 11. Let G be a graph such that A{G) has rank r, and let X C 
V{G) with \X\ = r. Then X is independent for A[G) iff dei G[X] = 1. 

Note that the independent sets X of cardinality equal to the rank are pre- 
cisely the bases of a matrix A. 

The following result is easy to see now from Proposition II II 

Lemma 12. Let G be a graph such that A{G) has rank r. Each element of J-q 
is of cardinality r. 

Proof. If there is an X € J-g of cardinality q > r, then the columns of ^(G[X]) 
are linearly independent, and thus so are the columns of A{G) corresponding to 
X. This contradicts the rank of A{G). 



1 Clearly, for a matrix A, detj4[X] 7^ implies that X is independent for A. The reverse 
implication is not valid in general. 
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Finally, assume that there is an X G Tq of cardinality q < r. Since the 
columns of A(G'[X]) are linearly independent, so are the columns of A{G) cor- 
responding to X. Since A{G) has rank r, X can be extended to a set X' 
with cardinality r. Hence by Proposition 1111 det G [X'] = 1 with X' Z) X — a 
contradiction of X G J-q- Q 

Example 13. We continue ExamplejH Let again G be the graph on the lower- 
left corner of Figure[3] Then the elements Fq — {{p, (?, s}, {p, 9, r}, {q, r, s}} are 
all of cardinality 3 — the rank of A{G). Moreover, Fa*{p,q] ~ {V} and \V\ — 4 
is equal to the rank of G * {p, q]. 

Combining Proposition II II and Lemma [T^ we have the following result. 

Corollary 14. Let G he a graph, and let X C V{G). Then X is a basis for 
AiG) tffXeTo. 

Equivalently, with Vg = {I, D) from Section[51 Corollary[T3]states that max(/) = 

By Corollaries [T] and [14] we have now the following. 

Lemma 15. Let G and G' be graphs. Then Tq = iff^cj{G) = kcr(G"). 

Recall that Lemma [7] shows that the dual pivot retains the kernel. We may 
now conclude from Lemma 1151 that also J-q is retained under dual pivot. It is 
the main result of this paper, and, as we will see in Section [3 has an important 
application. 

Theorem 16. Let G he a graph, and let X C V. Then — J-g*x if the 
right-hand side is defined. 

In particular, the case X — V, we have J-q+i = J'g-^+i if C! is invcrtible 
(over F2). 

Let Og = {GiX I X C V,det{G + I)[X] = 1} be the orbit of G under 
dual pivot, and note that G G Og- By Theorem [TCI if Gi,G2 G Og, then 
Jgi = J^G2- Note that the reverse implication does not hold: e.g. Oj = {/} 
and J-/ = {V}, while clearly there are many other graphs G with det G = 1 
(which means Tg — {V}). 

Example 17. We continue Example [S| Let again G be the graph on the lower- 
left corner of Figure [3| Then G*{p} is depicted on the lower-right corner of 
Figure[31 We have J^g*{p} = {{p, s}, {p, q, r}, {q, r, s}}, see Figure[ll so indeed 

— ^Gi{p}- 

For symmetric V x matrices A over F2, Theorem [TSI states that if A can 
be partially inverted w.r.t. Y (- V , where Y is maximal w.r.t. set inclusion, 
then this holds for every matrix obtained from A by dual pivot. 
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6. Maximal Contractions 

For a graph G, we define the contraction oi G on X C V with det G[X] = 1, 
denoted by G*\X, to be the graph {G*X)\X — the pivot on X followed by the 
removal of the vertices of X. Equivalently, contraction is the Schur complement 
applied to graphs. A contraction of G on X is maximal if there is no Y D X 
such that detG[y] = 1, hence if X £ J^q. The graph obtained by a maximal 
contraction on X is a discrete graph G' (without loops). Indeed, if G' were 
to have a loop e — {u} or an edge e — {u, v} between two non-loop vertices, 
then, since detG[X®e] = det((G * X)[e]) = det((G * \X)[e]) ^ 1, X ® e D X 
would be a contradiction of the maximality of X . Moreover, by Lemma [T^ the 
number of vertices of G' is equal to the nullity (dimension of the kernel, which 
equals the dimension of the matrix minus its rank) of G. 

Remark 18. In fact, it is known that any Schur complement in a matrix A has 
the same nullity as A itself — it is a consequence of the Guttman rank additivity 
formula, see, e.g., [l^. Section 6.0.1]. Therefore the the nullity is invariant under 
contraction in general (not only maximal contraction) . 

By Theorem 1161 we have the following. 

Corollary 19. The set of discrete graphs obtainable through contractions is 
equal for G and G*X for all X C V with det(G + I)[X] = I. 

In this sense, all the elements of the orbit Og have equal "behaviour" w.r.t. 
maximal contractions. 

Example 20. We continue the example. Recall that, from Example [T71 = 
■^G*{p} — {{Pi 'Zi {Pi '^}: {^i ''j s}}- The elementary contractions starting 
from G and G*{p} are given in Figure O Notice that the maximal contractions 
of G and G*{p} obtain the same set of (discrete) graphs. 

It is important to realize that while the maximal contractions (correspond- 
ing to Tq) are the same for graphs G and G*X, the whole set of contractions 
(corresponding to M.g) may be spectacularly different. Indeed, e.g., in Exam- 
pledOl the elementary pivots for G are *{q}, *{p,s}, *{p,r}, and *{r,s}, while 
the elementary pivots for G*{p} are *{r}, *{s}, and *{p,q} (see FigureH]). 

7. Application: Gene Assembly 

Gene assembly is a highly involved and parallel process occurring in one- 
cellular organisms called dilates. During gene assembly a nucleus, called mi- 
cronucleus (MIC), is transformed into another nucleus called macronucleus (MAC) 
Segments of the genes in the MAC occur in scrambled order in the MIC . Dur- 
ing gene assembly, recombination takes place to "sort" these gene segments in 
the MIC in the right orientation and order to obtain the MAC gene. The trans- 
formation of single genes from their MIC form to their MAC form is formally 
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Figure 5: Elementary contractions starting from G and G*{p}. 



modelled, see [1, [13, 0, as both a string based model and a (almost equivalent) 
graph based model. It is observed in [3| that two of the three operations in the 
graph based model are exactly the two elementary principal pivot transform 
(PPT, or simply pivot) operations on the corresponding adjacency matrices 
considered over F2. The third operation simply removes isolated vertices. 

Maximal contractions are especially important within the theory of gene 
assembly in ciliates — such a maximal sequence determines a complete trans- 
formation of the gene to its MAC form. We first recall the string rewriting 
system, and then recall the generalization to the graph rewriting system. 

Let A be an arbitrary finite alphabet. The set of letters in a string u over A 
is denoted by L{u). String u is called a double occurrence string if each x G L{u) 
occurs exactly twice in u. For example, u = 41215425 is a double occurrence 
string over L{u) = {1, . . . , 5}. Let A = {x \ x ^ A} with ACi A = 0, and let 
A = AlJ A. We use the "bar operator" to move from A to A and back from A 
to A. Hence, for x G A, x — x. For a string u — X1X2 ■ ■ ■ Xn with Xi € A, the 
inverse of u is the string u = XnXn-i ■ ■ ■ xi. 

We define the morphism || • || : (A)* — > A* as follows: for x G A, \\x\\ = x if 
X G A, and = x if x & A, i.e., ||a;|| is the "unbarred" variant of x. Hence, 
e.g., ||253|| = 253. A legal string is a string u e {A)* where ||m|| is a double 
occurrence string. We denote the empty string by A. 

Example 21. The string u — qpsqrpsr over A with A = {p,q,r,s} is a legal 
string. As another example, the legal string 3445675678932289 over B with 
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-8 = {2,3, ...,9} represents the micronuclear form of the gene cor resp onding to 
the actin protein in the stichotrichous cihate Sterkiella nova, see |l4i 0] . 



It is postulated that gene assembly is performed by three types of elementary 
recombination operations, called loop, hairpin, and double-loop recombination 
on DNA, see [15|. These three recombination operations have been modeled as 
three types of string rewriting rules operating on legal strings @, Q — together 
they form the string pointer reduction system. For all x,y G A with ||a;|| ||y|| 
we define: 

• the string negative rule for x by snrx{uixxu2) = uiU2, 

• the string positive rule for x by spr^{uixu2XU3) — U1U2U3, 

• the string double rule for x,y by sdrx.y{uixu2yu3XUiyu5) — U1U4M3W2U5, 

where ui,U2, • • ■ ,^5 are arbitrary (possibly empty) strings over A. 

Example 22. Let again u = qpsqrpsr be a legal string. We have spry{u) = 
sprpsr. And moreover, spr.p spr^ spr^(M) — ss. Finally, snr^ spr^ spr^ spr^(u) = 
A. 

We now define a graph for a legal string representing whether or not intervals 
within the legal string "overlap". Let u = X1X2 ■ ■ ■ Xn be a legal string with 
Xi & A for I < i < n. For letter y g let 1 < i < j < n be the 

positions of y in u, i.e., = — y. The y-interval of u, denoted by 

intvj^, is the substring XkXk+i ■ ■ ■ xi where k = i li Xi — y and fc = i + 1 if 
Xi = y, and similarly, I = j if Xj = y and I — j — 1 \i Xj = y (i.e., a border 
of the interval is included in case of y and excluded in case of y). Now the 
overlap graph of u, denoted by Qu, is the graph (V, -E) with V — L{\\u\\) and 
E = {{x,y} I X occurs exactly once in ||intvj,||}. Note that E is well defined 
as X occurring exactly once in the y-interval of u is equivalent to y occurring 
exactly once in the cc-interval of u. Note that we have a loop {x} G i? iff both 
X and X occur in u. The overlap graph as defined here is an extension of the 
usual definition of overlap graph (also called circle graph) from simple graphs 
(without loops) to graphs (where loops are allowed). See 113, Section 7.4] for a 
brief overview of (simple) overlap graphs. 

Example 23. The overlap graph Cy„ of m = qpsqrpsr is exactly the graph G of 
Example [3l 

It is shown in @,[l3], see also that the string rules snra;, spr^., and sdra;^^ 
on legal strings u can be simulated as graph rules gnr^, gpr^, and gdr^ ,^ on 
overlap graphs Gu in the sense that ^spr^(ti) — SP^xiSu), where the left-hand 
side is defined iff the right-hand side is defined, and similarly for gdr^ and 
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gnrjl. It was shown in [3| that gpr^ and gdr^ are exactly the two types 
of contractions of elementary pivots *\{x} and *\{x,y} on a loop {x} and an 
edge {x, y} without loops, respectively. The gnr^ rule is the removal of isolated 
vertex x. 

Example 24. The sequence spr^ spVp spr^ applicable to u given in Exam- 
ple [22] corresponds to a maximal contraction of graph G ~ Qu Example [3] as 
can be seen in Figure (5] 

Within the theory of gene assembly one is interested in maximal recombina- 
tion strategies of a gene. These strategies correspond to maximal contractions 
of a graph G (hence decomposable into a sequence ipi of contractions of elemen- 
tary pivots gpr and gdr applicable to (defined on) G) followed by a sequence 
ip2 of gnr rules, removing isolated vertices, until the empty graph is obtained. 
Here we call these sequences tp = Lp2'Pi of graph rules complete contractions. If 
we define the set of vertices v oi ip used in gnr^ rules by gnrdom((/3), then the 
following result holds by Corollarv ITOl 

Theorem 25. Let Gi,G2 G Og for some graph G, and let ip be a complete 
contraction of graph Gi . Then there is a complete contraction Lp' of G2 such 
that gnrdom(i^) = gnrdom(((5'). 

Hence, Theorem [25] shows that all the elements of Oq, for any graph G, have 
the same behaviour w.r.t. the applicability of the rule gnr^. 

A similar result as Theorem [25] was shown for the string rewriting model, 
see 0, Theorem 34jf|. It should be stressed however that Theorem [55] is real 
generalization of the result in [3] as not every graph has a string representation 
(i.e., not every graph is an overlap graph), and moreover it is obtained in a very 
different way: here the result is obtained using techniques from linear algebra. 

8. Discussion 

We introduced the concept of dual pivot and have shown that it has in- 
teresting properties: it has the same effect as the (regular) pivot and can be 
simulated by elementary row operations — consequently it keeps the kernel in- 
variant. The dual pivot in this way allows for an alternative definition of the 
(regular) pivot operation. Furthermore, we have shown that two graphs have 



■^There is an exception for gnr^.: although C?snr^(u) = S^^xiG^) holds if the left-hand side 
is defined, there are cases where the right-hand side is defined (x is an isolated vertex in Su) 
while the left-hand side is not defined {u does not have substring xx). This is why the string 
and graph models are "almost" equivalent. This difference in models is not relevant for our 
purposes. 

■^This result states that two legal strings equivalent modulo "dual" string rules have the 
same reduction graph (up to isomorphism). It then follows from [^ Theorem 44] that these 
strings have complete contractions with equal snrdom, the string equivalent of gnrdom. 
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equal kernel precisely when they have the same set of maximal pivots. From 
this it follows that the set of maximal pivots is invariant under dual pivot. 

This main result is motivated by the theory of gene assembly in ciliates in 
which maximal contractions correspond to complete transformations of a gene 
to its macronuclear form. However, as applying a maximal pivot corresponds to 
calculating a maximal partial inverse of the matrix, the result is also interesting 
from a purely theoretical point of view. 
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